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Abstract. Although there are now some tentative signs that the start of 
cycle 24 has begun there is still considerable interest in the somewhat unusual 
behaviour of the current solar minimum and the apparent delay in the true start 
of the next cycle. While this behaviour is easily tracked by observing the change 
in surface activity a question can also be asked about what is happening beneath 
the surface where the magnetic activity ultimately originates? In order to try to 
answer this question we can look at the behaviour of the frequencies of the Sun's 
natural seismic modes of oscillation - the p modes. These seismic frequencies 
also respond to changes in activity and are probes of conditions in the solar 
interior. 

The Birmingham Solar Oscillations Network (BiSON) has made measure- 
ments of low-degree (low-£) p mode frequencies over the last three solar cycles, 
and so is in a unique position to explore the current unusual and extended solar 
minimum. We compare the frequency shifts in the \ow-£ p-modes obtained from 
the BiSON data with the change in surface activity as measured by different 
proxies and show there are significant differences especially during the declining 
phase of solar cycle 23 and into the current minimum. We also observe quasi- 
biennial periodic behaviour in the p mode frequencies over the last 2 cycles that, 
unlike in the surface measurements, seems to be present at mid- and low-activity 
levels. Additionally we look at the frequency shifts of individual £ modes. 



1. Introduction 

The magnetic activity of the Sun is known to follow an approximately 11-year 
cycle ranging from quiet periods to much more active periods. Starting in March 
1755 these cycles have been labeled with a number and we are currently in the 
minimum between cycles 23 and 24. However, this current minimum is proving 
to be somewhat unusual, both in terms of how long it is lasting and just how 
quiet the Sun is, at least in terms of the more recently observed cycles. 

Surface measures of solar actively such as the number of visible sunspots 
and 10.7cm radio flux, all suggest that we are still in (or possibly just leaving) 
an extended minimum. Therefore, the question presents itself is there anything 
different about this minimum that we can learn from looking at the activity of 
the solar interior. 

This can be investigated by looking at the acoustic modes of oscillation of 
the Sun (known as p modes). Modes with low angular degrees (low-£) travel 
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deeply into the solar interior and so sample the conditions below the solar sur- 
face. Also, low-£ modes are truly global in nature and therefore are sensitive 
to a large fraction of the solar disc. It has been known for some time that the 
frequencies of these oscillations vary during the solar cycle, with the frequencies 
being highest when the magnetic activity is at a maximum. Therefore, by look- 
ing at the change in frequencies we can gain insight into cycle related processes 
that are occurring beneath the solar surface. 

The Birmingham Solar Oscillations Network (BiSON) has been collecting 
unresolved (Sun-as-a-star) Doppler velocity observations for over three decades. 
These types of observations are sensitive to the large scale 'global' oscillations 
of the sun, (i.e., those modes with the lowest angular degrees, £). Due to the 
fact the BiSON has observations going back into the 1970's, these data offer 
a unique opportunity to study the last three solar cycles. However, when the 
network was first established there was only one station and as such the early 
observations are sporadic and so the data quality is relatively poor. The quality 
and duty cycle of the data greatly improves after 1986 once three stations were 
operating and this allows us to study the last two solar cycles in their entirety. 



2. Method 



The precision with which one can measure the p-mode frequencies is directly 
related to the length of the time series. Therefore, one often wants to take 
as long a time series as possible in order to give the most precise estimates 
of the frequencies. However, when trying to track the change in frequencies 
over the solar cycle, one must strike a balance between length of time series 
and the time resolution of each frequency measurement. To this end an 8577- 
day BiSON time series starting on 14 April 1986 and ending on 7th October 
2009 was split into a set of 365-day time series. Each series was allowed to 
overlap the following by 91.25 days (i.e., a 4-time overlap), resulting in a total 
of 91 (non-independent) data sets. The power spectrum of each time series was 
parameterised using a modified Lorentzian model which w as fitted to the data 
using a standard likelihood maximisation method (eg., see lChaplin et al]|l999l ) 
enabling the mode frequencies to be determined. 

For each mode a minimum activity reference was taken from the fits to 
a data set lying on the boundary between cycle 22 and 23. The frequency 
shifts were then defined as the difference between this reference value and the 
frequency of the corresponding mode observed at different epochs. The size of 
the frequency shift has a known dependenc e on frequency and t his dependence 
was removed in the manner described in Chaplin et al.l ( 20071 ). A weighted 
average of the frequency shift for each time series was then determined, enabling 
the shift to be plotted over time. The average was made over modes ranging 
from 12 < n < 25 and < t < 3. 

A similar analysis was performed recently with BiSON data (see lBroomhall et al 



2009). Here we extend in time the results of that paper and additionally present 



the shifts of individual t modes. This allows us to make comments about the 
location of the solar variability, since each t has different latitudinal sensitivities. 



The New Solar Minimum. 
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Figure 1. Compariosn between the frequency shifts as determined from 
BiSON data and two different surface activity proxies. Left: The 10.7cm 
radio flux. Right: The International Sunspot number, both shown by the 
solid black line. 



3. Analysis and Results 

In this section we concentrate on the last two cycles (22 and 23) where the 
BiSON data are of sufficient quality to do meani ngful statistical tests. It is 
widely known (see Elsworth et al. 199Cll ; Howe 20081 . and references therein) that 



there is a good correlation between the shifts in the p-mode frequencies and 
other proxies that can be used to track the solar cycle. In Fig. [T] we compare 
the shifts with two of the more commonly used proxies, the 10.7cm radio flux 
(-^10. 7) an d the International Sunspot Number (ISN). The shifts show very good 
agreement with the 10.7cm radio flux up until around 2003, after which there is 
a significant disagreement relative to previous epochs. The agreement between 
the shifts and the ISN is not as good. This is thought to be because the ISN 
is predominantly sensitive to the strong component of the Sun's magnetic field, 
whereas the frequency shifts and radio flux are sensitive to both the strong and 
weak components. However, even with the poorer fit, it is still clear that the 
largest differences still occur over the latter part of cycle 23. 

Looking more closely at the final few points shown in Fig. [IJ it appears 
that there is a significant rise in the frequency shifts at the start of 2008, that 
may be pointing to the start of cycle 24. We do not see a similar rise in either 
the radio flux or the ISN. However, this rise in the shifts did not continue, and 
looking back at the boundary of the last cycle we can see there are similar rises 
and dips in the shifts throughout the last minimum too. 

If these dips and rises in the shifts are not associated with the main solar 
cycle then they do highlight an interesting notion in their own right. There are 
clearly bumps and dips at higher activity in both the shifts and other proxies, but 
if these continue into low activity periods too, then it might be a an indication 
of short-term pseudo periodic behaviour. 

To investigate this more thoroughly we have looked at the difference between 
heavily smoothed and unsmoothed data. The residuals of this comparison are 
shown in the left panel of Fig. [2) Short term periodic variations appear to 
occur on a scale of around two years. The variations seem to be stronger in 
the oscillations during the minimum periods but of similar magnitude to the 
radio flux during high activity. We can also look at a power spectrum of these 
residuals as shown in the right panel of Fig. [2j Although there is not a clear 
and obvious peak at short time periods there are a cluster of peaks around a 
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Residuals Power Spectrum 




Figure 2. Left: The residuals obtained when subtracting the unsmoothed 
frequency shifts (solid line) and radio flux (dotted line) from heavily smoothed 
values. Right: The power spectrum of these residuals. 



two to three year period. Similar "quasi-biennial" variation have been remarked 
upon before, s pecifically in the green co ronal emission line at 530. 3nm at high 
solar activity ( Vecchio &: Carbone 20081 ). Additionally, an explanation of such 
quasi-biennial behaviour has b een put forward by two different types of dynamos 
operating at different depths ( Benevolonskava lll998al lbl). 

We can gain further insight into what the shift in the frequencies are telling 
us by looking at the shifts of individual angular degrees. We show this for £ = 
- 3 in Fig. [3j The £ = shifts seem to fit the radio flux very well, with a reduced 
X 2 value of only 1.8 (which is considerably lower than the averaged £ value of 
3.5). This may be partly due to the fact that the £ = modes are more sensitive 
to the whole disc of the Sun than the higher £ modes. 

For the £ = 1 modes there is more structure than in the £ = case and this 
results in a poorer fit with the radio flux (the reduced x 2 is 3.3). The £ = 2 mode 
shifts show more similarity with £ = than £ = 1. However the fit to the radio 
flux is again poorer (a reduced x 2 01 3.1). This correlation between £ = and £ = 
2 modes and also £ = 1 and £ = 3 modes is an intriguing result. It might be due 
to the fact that correlated modes have somewhat similar latitudinal sensitivities. 
Both £ = and £ = 2 modes have an m = zonal component which is sensitive 
to higher latitudes, whereas £ = 1 and 3 modes only have sectoral components. 
Or alternatively, the effect could be an artifact of the fitting procedure since £ = 
and 2 modes and £ = 1 and 3 modes are fitted together in pairs. These ideas 
will be investigated further in the future. 

Again, we can look more closely at the final few points in each of the plots 
in order to investigate how the current solar minimum and start of cycle 24 is 
developing. For £ = we see a clear jump in the shifts at the beginning of 2008, 
but there is no obvious decline after this. For £ = 1 the rise at the start of 2008 
is much less pronounced than for £ = 0, and there is a definite decline after this. 
The £ = 2 shifts appear to start rising before the beginning of 2008 and continue 
upwards thereafter. Finally, for £ = 3 we see a short rise during 2008, although 
the error bars are sufficiently large that the rises and dips are not significant. 

Thes e results match quite closely with what is seen in an analysis of GOLF 
data (see Salabert et al.l l2009). However, in the GOLF data there was a clearer 
rise in the £ = and 2 modes and l ess of a dip in the £ = 1 data {£ = 3 modes 
were not included in this analysis). I Salabert et al. ( 20091 ) gave an explanation 
for this based on the latitudinal sensitivity of the modes. They suggest that, if 
it is assumed that at the time of the final few measurements the start of cycle 
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Figure 3. Compariosn between the frequency shifts as determined from 
BiSON data and the 10.7cm radio flux for each individual ^-mode in the 
range < I < 3 



24 had begun, then we would expect any increased activity to occur at higher 
latitudes first. The £ = and £ = 2 modes are more sensitive to these higher 
latitudes and hence may be seeing the effects of an increased activity due to the 
onset of cycle 24. The 1=1 modes, on the other hand, whose sensitivities are 
more confined towards the equator may not be seeing this increased activity yet. 

The main argument against this explanation is that the zonal m = com- 
ponents of the £ = 2 multiplet, which will be the component most sensitive to 
higher latitudes, only contributes about 10% to the determination of the fitted 
frequency. The frequency fit is actually much more dependent on the sectoral 
\m\ = 2 components, whose sensitivities are actually more concentrated around 
the equator than the £ = 1 sectoral components. Therefore it would seem that 
there is more investigation needed into the shifts of the individual £ modes to 
try and understand these differences. Also, as described above, we will inves- 
tigate the pair-by-pair nature of the fitting process to see whether it might be 
responsible for our observations. 



4. Summary and Discussion 

The shifts in oscillation frequencies over time is reasonably well correlated with 
other activity proxies such as the 10.7cm radio flux and ISN. However, there 
are some differences. The largest of which, at least over the two cycles we have 
data for, occur recently, starting on the downward part of cycle 23 and heading 
into the unusual and extended minimum between cycle 23 and 24. The current 
minimum in the frequency shifts is considerably deeper than in the proxies when 
compared with the previous minimum and there is clear structure in the shifts 
that does not appear in the proxies. Additionally, there is some evidence for 
pseudo-periodic short term (2-year) variations in activity on top of the 11-year 
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cycle. Frequency shifts show sharper amplitude variations during the quiet sun 
periods compared with the radio flux. 

Since the frequency shifts are sensitive to conditions below the surface it 
is likely that these differences are due to changes in the magnetic flux that 
are occurring in the solar interior. These changes then may have either yet to 
manifest on the surface, or are attenuated before ever reaching there. 

Alternatively the differences could be due to zonal effects, since the oscil- 
lations are not all equally sensitive at different latitudes of the Sun. We have 
looked at the shifts of modes with different I in order to gain a clearer picture of 
this. The I = modes show the closest match with the radio flux. The higher 
degree modes all show greater discrepancies. The strongest regions of magnetic 
flux tend to start at latitudes near the poles and then drift towards the equator 
as the cycle progresses. Therefore, modes that are only sensitive to ceratin lati- 
tudes will see their frequency shift's respond to this drift in magnetic flux. This 
will ultimately add further structure into the plots of the frequency shifts that 
may not be apparent in the I = modes or radio flux plots. 

If we assume that cycle 24 had started at the time of the last few points 
in the data set used here, then the shape of the plots for i = 0, 1, 2, and 3 
might also be explain ed by this latitudinal sensitivity argument as suggested by 
Salabert et al. ( 20091 ). Although the fact that the £ = 2 modes showed a similar 



upturn to the £ = modes in the final few observations, may actually work 
against this theory since the zonal m = component, which does have a greater 
sensitivity to higher latitudes, only has a small effect on the fitted frequency of 
the £ = 2 multiplet. 

It is clear that further work on the individual £ modes is needed to help 
clear up this matter. In the future we hope to continue the work on individual 
£ modes by comparing the frequency shifts with decomposed magnetic maps. 
This will allow us to compare the frequency shifts with the change in magnetic 
flux over time for similar zonal regions. 
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